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Introduction 


Purpose 


SWIFT/486  (Sandia  Waste-Isolation  Flow  and  Transport)  is  a  three- 
dimensional,  finite-difference  code  which  can  be  used  to  simulate  steady-state 
or  transient  flow  and  transport  of  chemicals  (including  brine  and  radionuclide) 
and  heat  in  porous  or  fractured  geologic  media.  The  geologic  media  may  be 
homogeneous,  isotropic,  heterogeneous,  and/or  anisotropic.  The  transport 
processes  which  may  be  modeled  by  SWIFT/486  include  advection,  disper¬ 
sion,  sorption,  decay,  and  leaching.  The  equations  for  fluid  flow,  heat  trans¬ 
port,  and  brine  transport  are  coupled  by  the  pore  fluid  velocity,  fluid  density, 
fluid  viscosity,  and  porosity.  Flow  and  mass  transport  in  fractured  media  are 
modeled  using  the  dual-porosity  approach  (fracture-matrix).  It  is  assumed  that 
the  transport  processes  in  rock-matrix  are  one  dimensional  in  a  lateral  direc¬ 
tion  relative  to  that  in  the  fracture  (Bear  and  Braester  1972,  Pruesses  and 
Narasimhan  1982).  Salt  dissolution  and  waste  leaching  algorithms  are  other 
optional  features  included  in  SWIFT/486.  The  purpose  of  the  waste  leach 
model  is  to  determine  the  source  rate  at  which  a  radionuclide  from  a  reposi¬ 
tory  is  dissolved  into  a  solution.  The  salt  dissolution  formulation  is  described 
in  detail  by  Nolen  et  al.  (1974).  Fluid  flow  of  variable  densities  and/or  vis¬ 
cosities  also  may  be  modeled  by  SWIFT/486.  Either  a  radial  or  Cartesian 
coordinate  system  can  be  used  for  domain  discretization.  The  present  version 
of  SWIFT/486  is  classified  as  a  single  phase  and  saturated  flow  model.  More 
detailed  information  on  the  code  capabilities  is  presented  by  Reeves  et  al. 
(1986a). 

The  SWIFT/486  program  is  an  enhanced  version  of  its  predecessor  codes, 
SWIP,  SWIPR,  SWIFT,  SWIFT  II,  and  SWIFT  III.  The  codes  SWIPR, 
SWIFT,  and  SWIFT  II  are  available  from  the  National  Technical  Information 
Service  and  the  Energy  Science  and  Technology  Center.  SWIFT  III  and 
SWIFT/486  are  available  from  GeoTrans,  Inc.,  Sterling,  Virginia.  Table  1 
provides  a  summary  of  available  references  on  each  of  the  above  versions  of 
SWIFT. 

The  PC  version,  SWIFT/486,  has  been  developed  for  the  Intel  80386  and 
Intel  80486  CPU  processors  using  the  FTN77/486  compiler  developed  by  the 
University  of  Salford.  An  optional  post  processing  program,  UNSWIFT, 
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Table  1 

Development  of  SWIFT/486 

Code 

Code  Developer 

Source  of 

Funding 

Reference 

SWIP 

Intercomp,  Inc. 

U.S.  Geological 
Survey  (USGS) 

Intercomp  (1976) 

SWIPR 

Intera,  Inc. 

USGS 

Intera  (1979) 

SWIFT 

Intera,  Inc. 

National  Reserch 
Council  (NRC) 

Oillon  et  al.  (1978),  Reeves 
and  Cranwell  (1981),  Finley 
(1981),  Ward  et  al.  (1984) 

SWIFT  II 

GeoTrans,  Inc. 

NRC 

Reeves  et  al.  (1986a), 

Reeves  et  al.  (1986b), 

Reeves  et  al.  (1986c) 

SWIFT  III 

GeoTrans,  Inc. 

GeoTrans 

Ward  (1987) 

SWIFT/386 

GeoTrans,  Inc. 

GeoTrans 

Ward  (1991) 

SWIFT/486 

GeoTrans,  Inc. 

GeoTrans 

Ward,  Harrover,  and  Vincent 
(1993) 

Source:  GeoTrans  (1993). 

provides  an  interface  to  the  contouring  program  SURFER  (Golden  Software, 
1993).  UNSWIFT  can  read  pressure,  temperature,  brine,  or  nuclide  concen¬ 
tration  data  from  a  SWIFT/486  output  file  (with  suffix  .MAP)  and  prepare  a 

grid  file  compatible  with  SURFER.  For  this  evaluation,  UNSWIFT  was  not 
tested. 


Scope  of  Report 

This  report  describes  the  evaluation  of  SWIFT/486  by  comparing  computed 
results  with  six  selected  analytical  solutions  for  several  flow  and  solute  trans¬ 
port  scenarios  of  varying  boundary  conditions  and  solute  sources  in  porous 
media.  The  analytical  solutions  were  selected  from  those  so-called  classic 
problems  such  as  Theis  (1935)  and  Hantush  (1960)  radial  problems  and  from 
the  latest  published  solutions  such  as  those  by  Batu  (1984)  and  Beljin  (1993). 
Some  of  the  solutions  were  also  given  in  the  SWIFT/486  documents.  The 
input  parameters  were  selected  from  the  SWIFT/486  reports  and  other  pub¬ 
lished  documents.  Tlie  analytical  solutions  used  are  useful  tools  to  test  and 
initially  verify  the  accuracy  of  SWIFT/486  algorithms  according  to  certain 
assumptions.  The  analytical  solutions  normally  are  exact  within  a  limited 
range  of  parameters.  These  solutions  are  easy  to  apply  and  require  fewer 
input  parameters.  Although  analytical  solutions  have  limited  applications  for 
real  field  problems,  they  have  been  used  extensively  to  check  the  correctness 
of  numerical  codes.  The  selected  analytical  solutions,  data  input,  and 
SWIFT/486  results  are  described  herein. 
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Model  Description 


Aquifer  Submodels 

Although  SWIFT/486  is  not  classified  as  a  multiphase  flow  and  transport 
code,  many  of  the  variables  in  the  SWIFT/486  code  and  terminology  used  in 
the  documentation  are  derived  from  petroleum  reservoir  engineering.  Hydro¬ 
geologists  unfamiliar  with  such  terminology  are  referred  to  Aziz  and  Settari 
(1979)  for  further  information.  Some  of  the  important  terms  used  in  SWIFT/ 
486  are  briefly  described  below.  Petroleum  reservoirs  are  usually  bounded 
partly  or  completely  by  the  water  saturated  zone  (aquifer).  In  SWIFT/486, 
the  simulation  region  is  divided  into  two  subregions,  an  inner  region  (reser¬ 
voir)  and  an  outer  region  (aquifer).  The  term  “reservoir”  is  applied  to  that 
portion  of  the  system  for  which  detailed  information  is  needed.  The  term 
“aquifer”  is  used  for  the  remaining  portion  of  the  system.  However,  hydroge¬ 
ologists  consider  both  the  inner  region  and  outer  region  as  aquifer  (saturated 
zone  of  water).  In  SWIFT/486,  the  aquifer  submodel  provides  boundary  con¬ 
ditions  for  the  reservoir.  A  diagram  showing  the  relationship  between  aquifer 
and  reservoir  for  a  cylindrical  geometry  is  shown  in  Figure  1  (Reeves  et  al. 
1986a).  The  inner  dimension  is  the  external  radius  of  the  reservoir.  The 
outer  dimension  is  the  external  radius  chosen  for  the  aquifer,  where  can 
have  a  finite  or  infinite  value.  The  default  value  for  in  SWIFT/486  is 
infinity.  The  aquifer  (outer  region)  thickness  A/z  is  input  into  SWIFT/486 
through  a  permeability-thickness  product  (average  value  of  transmissibility) 
and  a  porosity-thickness  product  defined  in  the  input  file.  Three  different 
aquifer  (boundary  conditions)  submodels  are  provided  in  the  SWIFT/486  pro¬ 
gram:  an  unsteady-state  aquifer,  a  steady-state  aquifer,  and  a  pot  (no-flow 
condition)  aquifer.  In  each,  a  type-three  condition  (e.g.,  Cauchy  boundary 
condition)  is  provided  for  each  boundary  grid  block.  The  rate  of  flow  from 
aquifer  to  reservoir  varies  with  the  pressure  change  within  the  block. 


Aquifer-influence  function  and  boundary  conditions 

Aquifer-influence  functions  are  analytic  submodels  used  for  treating  both 
external  and  internal  boundaries.  SWIFT/486  has  three  submodels  which  are 
coupled  to  the  reservoir  implicitly  under  the  influence  functions.  The  aquifer- 
influence  functions  are  designed  to  save  computational  time  by  simplifying 
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Figure  1 .  Geometric  parameters  of  the  reservoir  and  aquifer  in  cylindrical 
geometry  (after  Reeves  et  al.  1986a) 


simulation  in  the  peripheral  domain  regions  where  the  detailed  information  is 
not  needed.  Aquifer-influence  functions  can  be  used  to  specify  boundary  con¬ 
ditions  for  pressure,  temperature,  and  chemical  concentration.  The  default 
flow  boundary  condition  is  no  flow  across  the  aquifer-reservoir  boundaries. 


Initial  conditions 

SWIFTM86  assumes  initially  that  flow  in  a  hydrologic  system  is  in  hydro¬ 
static  equilibrium.  The  user  must  provide  an  initial  fluid  pressure  at  any 
chosen  reference  point  needed  to  calculate  hydrostatic  pressures  for  other 
poinp  in  the  domain  (Reeves  et  al.  1986a).  Constant  flux  of  flow  can  also  be 
specified.  The  initial  condition  for  temperature  is  specified  using  an  interpola¬ 
tion  function  (Reeves  et  al.  1986a).  The  initial  conditions  for  both  brine  and 
radionuclide  concentration  may  be  input  directly  by  the  user. 


Well  submodel 

A  source  or  sink  can  be  specified  in  SWIFT/486  using  the  well  option 
defined  in  the  code.  The  wells  also  may  be  used  to  simulate  both  a  nonpoint 
aquifer  recharge  from  upland  areas  and  an  aquifer  discharge  into  rivers  or 
streams.  This  model  option  can  be  used  to  switch  from  flow  rate  to  pressure 
control.  Well  option  is  also  used  to  specify  injection  (recharge)  or  pumping 
(discharge)  sources  of  fluid,  heat,  brine,  or  radionuclide. 

SWIFT/486  has  a  special  numerical  treatment  for  fluid  production/ 
injection  wells.  This  technique  has  been  used  in  petroleum  reservoir  engineer¬ 
ing.  For  each  well,  three  parameters  are  defined:  well  index,  mobility,  and 
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rate  allocation.  The  well  index  characterizes  the  transmitting  capability  of  the 
well  surrounding  a  region  (skin).  The  well  index  may  be  estimated  by  a  one¬ 
dimensional,  steady-state  flow  solution  which  is  a  function  of  the  hydraulic 
conductivity  of  well  surroundings.  The  horizontal  length  of  skin  may  be  cho¬ 
sen  as  the  size  of  the  block  (grid  cell).  In  this  case  the  well  index  is  related  to 
the  transmissibility.  The  mobility  parameter  is  similar  to  the  well  index 
except  it  is  a  function  of  thickness  of  aquifer  layers.  A  fractional  allocation 
factor  is  assigned  to  each  layer.  Each  factor  is  assumed  to  be  proportional  to 
the  thickness-permeability  product  for  that  layer.  The  rate  allocation  parame¬ 
ter  assigns  the  way  that  the  flow  partitioning  can  take  place  within  the  layers. 
The  rate  allocation  can  be  on  the  basis  of  mobilities  (flow  rates)  or  mobilities 
and  pressure  drops.  These  parameters  are  described  in  detail  by  Reeves  et  al. 
(1986a,b,  and  c). 


Global,  Local,  and  Primary  Equations 

SWIFT/486  may  be  used  to  simulate  flow  and  mass  transport  in  both 
porous  and/or  fractured  media.  Two  separate  sets  of  equations  are  used:  one 
for  porous  or  fractured  zones  and  a  second  set  to  describe  flow  and  transport 
in  matrix  blocks  of  a  fractured  zone.  The  equations  which  are  used  for 
porous  media  or  fractured  zone  are  called  “global”  equations.  The  equations 
which  describe  flow  and  mass  transport  in  the  matrix  part  of  a  fractured  zone 
are  called  “local”  equations.  The  term  “global”  or  the  term  “local”  is  used  in 
SWIFT/486  for  the  parameters  calculated  by  the  related  equations.  In 
SWIFT/486,  primary  equations  refer  to  equations  for  flow,  heat,  and  brine 
transports  because  density,  viscosity,  porosity,  and  enthalpy  are  functions  of 
pressure,  temperature,  and  brine  concentration.  Two  steady-state  solutions 
options  are  provided  for  the  global  flow  and  brine  transport  equations.  Heat 
transport  and  radionuclide  transport  are  not  included  in  the  steady-state  option. 
The  code  will  permit  a  steady-state  solution  of  the  primary  equations  (flow  or 
brine)  with  the  transient  solution  of  radionuclide  or  heat  transport. 


Mass  Balance  Calculation 

SWIFT/486  calculates  local  and  global  mass  balances  for  flow  and  solute 
transport  including  heat,  brine,  unleached  radionuclide,  leached  but  not  dis¬ 
solved  radionuclide  (enhanced  transport),  dissolved  radionuclide  and  radio¬ 
nuclide  in  the  matrix  subsystem.  The  control  parameter  for  mass  balance 
output  is  named  LMBAL  on  record  M-2  in  the  input  file.  In  addition,  the 
model  user  can  control  the  mass  balance  output  for  both  the  global  system  and 
the  local  subsystems  by  specifying  the  parameter,  101,  in  the  R2-13  block  of 
the  input  file.  Mass  balance  at  any  time  of  a  simulation  also  can  be  written 
into  a  file  assigned  for  Unit  17. 
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Numerical  Considerations 


The  numerical  results  of  SWIFT/486  mass  transport  simulations  are  usually 
a  function  of  two  dimensionless  numbers,  Peclet  number  and  Courant  number. 
The  calculated  values  of  these  two  numbers  will  appear  on  the  screen  during  a 
simulation  and  in  the  main  output  file.  Because  of  the  significance  of  Peclet 
and  Courant  numbers  and  their  effects  on  the  SWIFT/486  simulation  results, 
these  numbers  are  defined  here.  The  Peclet  number  is  the  ratio  of  advective 
to  dispersive  transport,  and  in  the  x-direction  is  given  as  follows: 

p  =  ^  ^  ^ 

'  ^  a  a 


where 

=  velocity  in  x-direction  [LT"’] 

a  =  longitudinal  dispersivity  [L] 

Av  =  grid  spacing  in  x-direction  [L  ] 

=  dispersion  coefficient  in  x-direction  [L'T"'] 

The  Courant  number  is  the  ratio  of  a  distance  travelled  by  constituent  (e.g, 
contaminat,  solute,  etc.)  within  a  time  step  A/  to  the  grid  dimension  in  a  flow 
direction.  For  flow  in  the  x-direction,  the  courant  number  Co  is  defined  as 
follows:  "" 


The  Peclet  number  and  Courant  number  are  functions  of  grid  spacing  time 
step,  and  flow  velocity.  In  SWIFT/486,  the  block-size  and  time-step  restric¬ 
tions  are  not  overly  severe  for  many  problems  (Reeves  et  al.  1986a).  In  some 
cases,  the  convection  terms  (velocity  terms)  in  the  transport  equations  may 
cause  some  numerical  errors.  These  errors  are  introduced  into  the  solution  as 
numerical  dispersion  and  the  overshoot-undershoot  phenomena.  To  reduce 
these  numerical  errors,  the  convection  terms  in  transport  equations  usually  are 
modified.  Several  techniques  have  been  developed  to  overcome  these  pro¬ 
blems.  Among  these  are  the  method  of  characteristics  (Carder  et  al.  1964, 
Bredehoeft  and  Binder  1973),  a  higher  order  Galerkin  method  (Price  et  al. 
1968,  Binder  1973],  various  upstream-weighting  and  asymmetric-weighting 
strategies  (Nolen  and  Berry  1972,  Christie  et  al.  1976),  and  the  distributed- 
velocity  method  (Campbell  et  al.  1981). 
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The  SWIFT/486  code  has  several  options  for  treating  numerical  dispersion 
and/or  overshoot-undershoot  problems.  In  Table  2,  v  is  a  generalized  Darcy 
flow  velocity;  Ax,  Ay,  and  At  are  grid  spacing  in  x-direction,  y-direction,  and 
time,  respectively;  is  the  longitudinal  dispersivity;  u  is  the  Darcy  velocity; 

is  molecular  difftision  coefficient;  4>  is  the  porosity;  is  the  formation 
density;  is  the  adsorption  coefficient;  is  the  heat  conductivity  parame¬ 
ter;  Cp  is  the  specific  heat  of  the  fluid;  is  the  specific  heat  of  the  rock 
(formation);  and  is  the  compressibility  of  fluid. 


Table  2 

Numerical  Criteria  for  Brine,  Heat,  and  Radionuclide  Transport 
(Reeves  et  al.  1986a) 

Scheme^ 

Numerical  Dispersion 

Dispersion  Criterion 

Overshoot  Criteria 

CIT-CIS 

None 

None 

1^// Aa'+2Z7A//Aa^  <  2 
vt^x!  2<  D 

CIT-BIS 

vtiXl2 

vhxl  2«  D 

vtt! Lx^2DLtl <2 

BIT-CIS 

v^Ltl2<D 

vhx!  2<D 

BIT-BIS 

vLxl2*v^htl2 

vhxl2*v^ttl2<D 

None 

D=  !  (K<l>)\  ,  radionuclide  or  brine  transport 

y  <7  heat  transport 

'  CIT  =  centered  in  time;  BIT  =  backward  in  time;  CIS  =  central  in  space;  BIS  =  backward  in  space^ 

Input  File 

Two  important  elements  of  an  input  file  in  SWIFT/486  consist  of  data  on 
geometric  gridding  and  on  numerical  criteria  parameters.  In  regard  to  geo¬ 
metric  gridding,  both  three-dimensional  Cartesian  (x,y,z)  and  axially  sym¬ 
metric  coordinates  (r,z)  can  be  modeled  by  SWIFT/486.  Discretization  in  a 
Cartesian  system  is  done  through  direct  input  of  the  increments.  The  user 
may  generate  a  mesh  system  by  specifying  all  values  of  increment  in  the  x-, 
y-,  and  z-direction  (DX,  DY,  and  DZ)  in  the  input  file.  For  the  radial  (r,z) 
coordinate,  the  mesh  may  be  generated  automatically  by  assigning  a  few 
parameters;  alternatively,  all  mesh  data  may  be  defined  explicitly  by  the  user. 

The  discrete  geometry,  in  either  the  Cartesian  coordinate  system  or  cylin¬ 
drical  coordinate  system,  is  called  a  global  in  order  to  distinguish  it  from  local 
discretization.  For  fractured  zone  discretization,  some  grid  blocks  may  be 
defined  as  dual-porosity  or  doubly  porous  media.  For  such  blocks,  a  local 
mesh  is  automatically  generated.  The  numerical  criteria  are  controlled  by 
parameters  assigned  in  the  input  file  by  the  user.  These  criteria  control 
numerical  dispersion,  overshot-undershot  errors,  adjustment  of  decay 
constants,  etc.  Detail  on  forming  an  input  file  is  given  in  Ward  et  al.  (1993). 
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Output  File 


The  main  output  file  created  by  running  SWIFT/486  has  the  suffix  .OUT. 
Other  auxiliary  files  (Table  3)  may  also  be  generated  for  other  purposes  by 
assigning  several  output  control  parameters  in  the  input  file.  Table  3  shows  a 
list  of  available  output  files  in  SWIFT/486.  The  main  output  file  and  its  unit 
number  are  shaded  in  Table  3.  For  more  detail,  the  reader  is  referred  to 
SWIPT/486  user’s  manual  (Ward  et  al.  1993). 


Table  3 

Available  Input/Output  File  Options 

Unit 

Function 

Default  File  Suffix 

4 

Input  for  restart  calculation 

.RST 

7 

Output  for  streamline  postprocessing 

.VL 

8 

Output  for  subsequent  restart  calculation 

.WR 

9 

Output  for  nuclide  monitor  post  processing 

.NM 

10 

Output  for  contouring  based  on  mapping  options 

.XYZ 

11 

Input  for  heterogeneous  reservoir  R1-21 

.BIN 

12 

Input  and  output  for  plotting  via  SWIFT 

.WL 

13 

Output  for  contour  mapping  using  MODFLOW  format 
lUNSWIFT  program  reads  this  file) 

.MAP 

15 

Standard  80  column  input 

.DAT 

■lej. 

-Standard  132  column  output  for  printer  . :: 

.OUT 

17 

Output  for  mass  balance  summary 

.MBL 

18 

Output  for  an  aquifer  influence  function  flux  values 

.AIF 
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A  numerical  code  such  as  SWIFT/486  may  be  evaluated  for  efficiency  of 
coding  by  checking  its  speed  of  running  (CPU  time)  for  specific  computer 
type  and  simulating  problems,  optimal  use  of  computer  storage,  convenience 
of  input/output,  program  portability,  and  diagnostic  messages.  In  this  investi¬ 
gation,  SWIFT/486  was  evaluated  for  all  the  above  evaluation  steps  except 
optimal  use  of  computer  storage.  The  numerical  accuracy  of  SWIFT/486  was 
evaluated  by  a  comparison  of  simulation  results  with  analytical  solutions.  The 
selected  analytical  problems  include  a  variety  of  initial  and  boundary  condi¬ 
tions.  The  model  was  reviewed  for  efficiency  of  coding  by  checking  its  speed 
of  run,  convenience  of  input/output  by  checking  data  input  and  output  infor¬ 
mation,  program  portability  by  using  it  in  two  different  computer  systems 
(DOS  and  UNIX),  and  available  diagnostic  messages  by  observing  those 
received  during  the  simulations.  Although  emphasis  in  this  report  is  given  to 
model  performance  against  the  analytical  solutions,  the  concluding  remarks  on 
all  of  the  above  evaluation  steps  are  provided  for  the  reader. 


Platform  for  Evaluation 

SWIFT/486-Version  2.53  has  been  designed  for  a  PC,  specifically  the 
80386  and  80486  processors,  offering  a  run-time  monitor  to  facilitate  the  pro¬ 
gress  and  status  of  batch  processing.  It  requires  the  FTN77/486  (FTN77/x86) 
Fortran  compiler,  version  2.6  and  higher.  The  simulations  described  here 
were  performed  on  a  486  PC  running  at  66  MHz  using  DOS  with  the  FTN77/ 
x86  Fortran  compiler.  The  use  of  this  specific  compiler  does  not  mean  an 
endorsement  by  the  U.S.  Army  Engineer  Waterways  Experiment  Station 
(USAEWES)  nor  by  any  other  branch  of  the  U.S.  Government. 


Modifications  for  the  Evaluation 

SWIFT/486  reads  input  parameters  from  a  formatted  file.  For  large  data 
input,  this  could  be  a  cumbersome  task.  To  provide  an  easy  way  to  enter  data 
in  the  input  file,  the  READ  statements  in  the  SWIFT/486  source  files  were 
modified  to  read  unformatted  parameters.  The  source  programs  then  were 
complied  and  linked  using  FTN/x86  supported  by  Salford  Software  (1993). 
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Example  Problems 


In  this  section  the  applications  of  SWIFT/486  to  six  problems  for  several 
flow  and  mass  transport  scenarios  are  discussed.  To  check  the  accuracy  of 
SWIFT/486  algorithms  within  a  limited  range  of  available  data,  the  simulated 
results  are  compared  with  the  analytical  solutions  for  simplified  groundwater 
problems.  Table  4  shows  a  summary  of  test  scenarios. 


Table  4 

Summary  of  Test  Scenarios 

Problem 

Number 

Description 

Coordinate 

System 

Reference 

1 

Fully  penetrating  well 

Radial 

Theis  (1935) 

2 

Fully  penetrating  well  in  a  leaky  aqui^ 
fer 

Radial 

Hantush  {I960.  1961) 

3 

Steady-state  horizontal  flow  in  a 
heterogeneous  aquifer 

Cartesian 

Batu  (1984) 

4 

Transport  in  a  plane  flow 

Radial 

Beijin  (1991,  1993) 

5 

Transport  from  a  continuous  source 

Cartesian 

Beijin  (1991.  1993) 

6 

Transport  of  a  solute  slug 

Cartesian 

Beijin  (1991.1993) 

Problem  1,  Fully  penetrating  well  in  a  confined  aquifer 

Objectives.  The  purpose  of  this  exercise  is  to  test  the  pressure  solution, 
aquifer-influence  function,  radial  geometry,  and  the  rate-controlled  well 
parameter. 

Problem  statement.  A  fully  penetrating  well  tapping  an  infinite,  homoge¬ 
neous,  confined  aquifer  is  pumped  at  a  constant  rate.  The  resulting  drawdown 
can  be  calculated  using  an  analytical  solution  developed  by  Theis  (1935).  The 
aquifer  storativity  S  and  transmissivity  T  are  calculated  using  rock  compress¬ 
ibility  c^,  the  water  compressibility  the  hydraulic  conductivity  K,  and  the 
aquifer  thickness  b  as; 

S  =  pg(f>  (c^  +  cjb 


where  g  is  the  gravitational  constant  (9.806  m/s)  and  other  parameters  are 
defined  in  Table  5.  The  water  compressibility  is  assumed  to  be  zero  =  0). 

Input.  The  specific  input  parameters  for  this  problem  are  given  in 
Table  5.  These  data  are  referenced  by  Ward  et  al.  (1984)  and  Ross  et  al. 
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Table  5 

Input  Parameters  for  Example  Problem  1  (Theis  Problem) 

Parameter 

Symbol/Unit 

Value 

Storativity 

5  (dimensionless) 

0.001 

Transmissivity 

T  (m*/s) 

0.001 

Pumping  rate 

Q  (m^/s) 

0.003 

Porosity 

<p  (dimensionless) 

0.2 

Hydraulic  conductivity 

K  (m/s) 

3.28  X  10"* 

Viscosity 

(Pa-s) 

0.001 

Density 

p  (kg/m^  ) 

999.5 

Rock  compressibility 

c,  (1/Pa) 

1.67  X  10  ' 

Aquifer  thickness 

b  (m) 

3.05 

Wellbore  radius 

(m) 

0.1143 

Aquifer  radius 

r.  (m) 

6096 

(1982)  and  are  used  here  in  both  SWIFT/486  simulations  and  the  analytical 
solutions. 


Numerical  specification.  The  SWIFT/486  radial  (cylindrical  geometry) 
coordinate  system  option  was  used  by  specifying  50  elements  in  the  radial 
direction  and  one  element  in  the  z-  and  y-directions.  The  well  radius  was 
=  0.1 143  m,  and  the  radius  of  the  domain  was  r,  =  6096  m.  The  selected 
domain  radius  satisfies  the  condition  of  infinite  domain  required  by  the  analyt¬ 
ical  solution.  The  discretization  in  the  radial  direction  was  done  automatically 
by  the  SWIFT/486  program.  (See  the  input  and  the  output  for  this  problem  in 
Appendix  A.) 

Output.  Output  for  this  example  (Figures  2  and  3  and  Appendix  A) 
includes  pressure,  time,  and  distance.  Because  the  analytical  results  were  in 
terms  of  drawdown  and  not  in  pressure,  the  output  from  SWIFT/486  was  con¬ 
verted  into  the  drawdown  s  [L]  by  the  following  equation: 


Pb-P  (r.  z,  t) 

P8 


(4) 


where 
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Py  =  specified  pressure  value  at  the  boundary.  Pa 

p(r,z,t)  =  calculated  pressure  at  a  radial  distance  r,  vertical  distance  z,  and 
time  t 

p  =  water  density 

g  =  gravitational  constant  (9.806  m/s) 

Results.  Graphical  comparison  of  numerical  approximation  and  analytical 
solution  are  presented  in  Figure  2  as  the  drawdown  versus  time  at  a  radial 
distance,  r  =  100  m,  and  drawdown  versus  distance  at  time,  t  =  100  days, 
respectively.  There  is  good  agreement  between  the  numerical  results  and  the 
analytical  results  at  all  distances  and  times. 

This  test  shows  that  SWlFT/486  correctly  calculates  spatial  and  temporal 
pressure  variations  in  a  radial  (cylindrical)  coordinate  system.  The  test  also 
verifies  that  the  aquifer-influence  function  works  according  to  its  purpose 
which  is  to  provide  boundary  conditions  for  an  infinite  domain. 


Problem  2,  Fully  penetrating  well  in  a  leaky  aquifer  with  storage 

Objectives.  The  objectives  in  Example  Problem  2  are  to  test  the  pressure 
solution,  the  coupling  of  vertical  flow  in  an  aquitard  with  horizontal  flow  in 
an  aquifer,  a  rate-controlled  well  condition,  and  the  aquifer-influence  function 
in  a  radial  coordinate  system. 

Problem  statement.  A  fully  penetrating  well  in  an  infinite,  homogeneous, 
and  isotropic  aquifer  is  pumped  at  a  constant  rate.  The  aquifer  is  bounded 
below  by  an  impervious  layer  and  above  by  a  semi-impervious  layer  or  aqui¬ 
tard.  Initially,  the  water  pressure  is  uniform  in  the  aquifer  and  the  aquitard. 
Flow  is  predominantly  vertical  in  the  aquitard  and  horizontal  in  the  aquifer. 
The  analytical  solution  for  thjs  problem,  assuming  homogeneous  and  isotropic 
aquifer  and  constant  pumping  rate,  was  given  by  Hantush  (1960,  1961). 

Input.  The  major  input  parameters  used  in  this  problem  (Table  6)  are 
described  by  Ross  et  al.  (1982)  and  Ward  et  al.  (1984). 

Numerical  specification.  The  SWIFT/486  radial  coordinate  system  option 
(cylindrical  geometry)  was  used  by  specifying  50  elements  in  the  radial  direc¬ 
tion,  one  element  in  the  y-direction,  and  two  elements  in  the  z-direction.  The 
well  radius  is  0.1143  m  and  the  radius  of  domain  r,  is  6,096  m.  The  selected 
domain  radius  satisfies  the  condition  of  infinite  domain  imposed  by  the  analyt¬ 
ical  solution.  The  discretization  in  the  radial  direction  was  done  automatically 
by  the  SWIFT/486  program  by  assigning  RI  =  0.2957  m  in  the  Rl-22  record 
in  the  input  file  (Appendix  A). 
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a)  Fixed  radial  distance  and  variable  time  (results  at  100  m  from 
the  well) 


0 

SWIFT 

THEIS 

J - - - L 


0  1 000  2000  3000  4000  5000  6000 

Distance  (m) 

b)  Fixed  time  and  variable  radial  distance  (results  at  time  = 

1 00  days) 


Figure  2.  Comparison  of  SWIFT/486  simulation  results  with  analytical  solu¬ 
tions  for  Theis  radial  problem 


Output.  The  SWIFT/486  output  is  given  in  Appendix  A.  From  that  out¬ 
put,  the  pressure  data  were  converted  into  the  drawdown  and  were  plotted 
versus  radial  distance  in  Figure  3. 

Results.  The  results  of  the  SWIFT/486  simulation  and  analytical  solution 
at  a  radial  distance  of  20  m  are  shown  in  Figure  3.  A  comparison  between 
these  results  indicates  that  there  is  a  good  agreement  between  the  numerical 
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Table  6 

Input  Parameters  for  Example  Problem  2  (Hantush  Problem) 

Parameter 

Symbol/Unit 

Value 

Aquifer  storativity 

S  (dimensionless) 

10“ 

Aquifer  transmissivity 

T  (m^/s) 

10^ 

Aquitard  specific  storativity 

S  (dimensionless) 

3.0  X  10^ 

Aquitard  hydraulic  conductivity 

K'  (m/s) 

3.0  X  10  '“ 

Aquitard  thickness 

b'  (m) 

Pumping  rate 

Q  (m^/s) 

0.014 

Aquitard  porosity 

0' 

0.4 

Aquifer  porosity 

'  0 

0.004 

Fluid  viscosity 

fJ  (Pa-s) 

0.001 

Fluid  density 

p  (kg/m^) 

1,000 

Rock  compressibility 

C/j(Pa  ^) 

7.67  X  10  ' 

Wellbore  radius 

0.1143 

Aquifer  radius 

^(m) 

6096.0 

Figure  3.  Results  of  SWIFT/486  simulation  and  leaky  aquifer  analytical 
solution  |R  =  20  m) 
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and  analytical  solutions  up  to  about  10  min  after  the  pumping  started.  The 
results  after  10  min,  however,  indicate  that  SWIFT/486  overestimated  the 
drawdown  compared  with  the  analytical  solution.  This  overestimation  of  the 
results  by  SWIFT/486  is  due  to  an  increase  of  the  time  step  as  the  simulation 
proceeded  with  time.  Similar  agreement  also  was  reported  by  Ward  et  al. 
(1984). 

This  test  verifies  the  accuracy  of  the  pressure  solution  for  a  leaky  aquifer 
with  radial  geometry  using  the  rate-controlled  well  condition  and  aquifer- 
influence  function  options. 


Problem  3,  Steady-state  horizontal  flow  in  a  heterogeneous  aquifer 

Objectives.  The  purpose  of  this  simulation  is  to  test  the  capability  of 
SWIFT/486  for  simulating  flow  through  a  simple  heterogeneous  porous 
medium  in  a  Cartesian  coordinate  system. 

Problem  statement.  A  steady-state,  one-dimensional,  horizontal  flow 
occurs  in  a  confined  and  rectangular  aquifer.  The  horizontal  hydraulic  con¬ 
ductivity  varies  continuously  as  an  exponential  function  of  x  and  z  (Batu 
1984).  The  vertical  hydraulic  conductivity  is  assumed  to  be  zero  everywhere. 


The  governing  equation  is  given  by: 


d  KJ^x,z) 
d  X 


d  h 
d  X 


+  K^ix,z) 


^  h 
a  ;r- 


=  0 


(5) 


in  which  the  hydraulic  conductivities  in  x-direction  and  z-direction  are 
defined  as: 


K^(x^)  =  a  exp  {bx  +  cz  +  d) 


K,(pc^)  =  0 


(6) 


where  a  [LT%  b  [L  'J,  c  [L"'],  and  d  [L\  are  arbitrary  constants. 

The  solution  for  Equation  5  (Batu  1984)  for  a  simulated  domain  of  hori¬ 
zontal  length  dimension  L  and  constant  boundary  heads  at  the  upstream  and 
the  downstream  /i2  is: 

h  Qc)  =  [exp  (  -  bL)  -  1]“'  [/i,  exp  (  -  bL)  ^ 

-  /j,  +  (/i,  -  exp  (  -  bx)\ 


Using  the  Darcy  flow  equation,  the  seepage  velocity  component  in  the 
x-direction  is  given  as: 
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u 


(8) 


0 


ah  [1  -  exp  (  -  bL)]'’(hj  -  exp  (cz  +  d) 


where  (f)  is  the  porosity. 

Input.  The  input  parameters  for  this  problem  were  taken  from  Batu 
(1984)  and  are  listed  in  Table  7. 


Table  7 

Input  Parameters  for  Example  Problem  3  (Batu  Problem) 

Parameter 

Value 

a 

0.006  m/s 

b 

0.05  m  ' 

c 

0.20  m  ’ 

d 

-2 

Porosity,  n 

0.4 

Horizontal  length,  L 

20  m 

Upstream  piezometric  head, 

11m 

Downstream  piezometric  head, 

6  m 

Numerical  specification.  The  domain  for  this  problem  has  a  dimension  of 
20  m  in  the  x-direction,  0.5  m  in  the  y-direction,  and  0.5  m  in  the  z-direction. 
A  total  of  40  elements  in  the  x-direction,  1  element  in  the  y-direction,  and 
10  elements  m  the  z-direction  were  used  for  this  application.  The  boundary 
conditions  consist  of  fixed  total  hydraulic  heads  at  x  =  0  and  x  =  20  m  and 
zero  fluxes  at  the  upper  and  lower  parts  of  the  domain  (Figure  4).  The  ’ 
hydraulic  conductivity  varies  in  the  x-  and  z-direction;  a  contour  plot  of  the 
hydraulic  conductivity  is  given  in  Figure  5. 

Output.  Output  for  Example  Problem  3  is  given  in  Appendix  A.  Seepage 
velocity  and  hydraulic  conductivity  data  are  plotted  in  Figure  6. 

Results.  As  shown  in  Figure  6,  the  numerical  results  are  almost  identical 
to  the  analytical  solution  results.  This  test  showed  that  SWIFT/486  is  capable 
of  solving  a  simplified  heterogeneous  porous  media  problem. 
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Figure  4.  Schematic  of  domain  with  boundary  conditions  specified  along 
the  domain 


Problem  4,  Transport  in  a  plane  radial  flow  field 

Objectives.  The  main  objective  in  Example  Problem  4  is  to  evaluate  the 
brine  transport  option  of  SWIFT/486  by  applying  it  to  solute  transport  of  a 
chemical  with  constant  water  density.  Other  objectives  are  to  check  the  ability 
of  SWIFT/486  to  calculate  the  velocity  field  around  an  injection  well,  to  simu¬ 
late  solute  transport  in  radial  flow,  and  to  use  the  aquifer  influence  function. 

Problem  statement.  This  problem  describes  the  dispersion  of  a  conserva¬ 
tive  solute  injected  via  a  fully  penetrating  well  in  a  confined,  homogeneous 
and  isotropic  aquifer.  The  main  assumptions  are  that  (a)  the  injection  rate  of 
fluid  is  constant,  (b)  the  regional  groundwater  velocity  is  negligible  compared 
with  the  velocity  created  by  injection,  and  (c)  steady-state  flow  occurs.  Note 
that  for  steady-state  plane  radial  flow,  the  product  of  velocity  times  radial  dis¬ 
tance  Vr  remains  constant. 

The  governing  advective-dispersive  equation  is  given  as  (Hoopes  and 
Harleman  1967): 

I  Alnr^]  - 

r  dr  \  dr)  dr 

where 

r  =  radial  distance  [L] 

D  =  dispersion  coefficient  [L-  T'] 

C  =  concentration  [M  L’^J 

V  =  average  flow  velocity  [L  T  'J 


dC 

dt 


(9) 
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Figure  5.  Contours  of  hydraulic  conductivity  as  function  of  x  and  z  (m/s) 


Figure  6.  Simulated  SWIFT  results  and  analytical  solutions  for  inhomoge¬ 
neous  domain  (Batu  1984) 


Input.  The  main  input  parameters  and  analytical  solution  for  this  problem 
are  from  Beijin  (1991,  1993).  A  summary  of  data  used  in  SWIFT/486  and 
the  analytical  solution  included  in  SOLUTE  software  package  by  Beijin  (1993) 
is  given  in  Table  8. 

Numerical  speciHcation.  A  cylindrical  domain  with  100  elements  in  the 
radial  direction  and  1  element  in  the  z-direction  was  specified  for  this  prob¬ 
lem.  The  radial  distance  was  50  m,  and  the  total  depth  was  3.048  m.  The 
radius  of  the  well  was  0. 1 143  m.  The  brine  transport  option  of  SWIFT/486 
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Table  8 

Input  Parameters  for  Example  Problem  4  (Beijin  Radial  Problem) 

Parameter 

Value 

Well  recharge  rate  Q 

25.0  m^/d  (4.59  gpm) 

Thickness  of  aquifer,  b 

3.048  m  (10.0  ft) 

Porosity,  n 

0.25 

Lateral  dispersivity,  Oj 

The  same  as  longitudinal 

Longitudinal  dispersivity, 

0.300  m  (0.984  ft) 

0.150  m  (0,492  ft) 

0.015  m  (0.049  ft) 

Time,  t 

20.0  d 

was  used  to  simulate  a  conservative  solute  transport  for  three  dispersion 
coefficients. 

Output.  Output  for  Example  Problem  4  is  included  in  Appendix  A.  The 
graphic  representations  of  scaled  or  normalized  (dimensionless)  concentration 
versus  the  radial  distances  at  a  time  of  20  days  for  different  dispersion  coeffi¬ 
cients  are  plotted  in  Figures  7-9. 

Results.  SWIFT/486  was  simulated  for  three  different  dispersivity  coeffi¬ 
cients  of  0.3  m,  0.15  m,  and  0.015  m.  As  shown  in  Figures  7  and  8,  there  is 
a  good  agreement  between  the  simulated  results  and  analytical  solutions  for 
the  two  former  dispersivities.  For  dispersivity  of  0.015  m  (higher  Peclet 
number),  the  simulated  results  deviated  slightly  from  the  analytical  solution  at 
certain  locations  (Figure  9). 

For  radial  problems,  SWIFT/486  incorrectly  calculates  the  Peclet  number.’ 
Therefore,  for  the  radial  problems,  the  user  can  ignore  a  warning  sign  of  a 
high  Peclet  number  that  may  appear  on  screen  or  in  the  output  file.  In  other 
words,  for  a  radial  problem,  even  if  the  Peclet  number  is  in  the  acceptable 
range,  the  computer  program  will  still  warn  the  user  of  a  high  Peclet  number. 

This  test  verified  that  the  “brine”  transport  options  of  SWIFT/486  worked 
correctly  to  solve  well  injection  of  a  dissolved  chemical  for  a  simple  problem. 


’  Personal  Communication,  (1994).  D.  Ward,  Vice  President,  GeoTrans,  Inc.,  Sterling,  VA  20166. 
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Figure  7.  Simulated  SWIFT/486  results  and  analytical  solutions  for  a  radial 
solute  transport  problem  with  dispersivity  coefficient  of  0.3  m 


Figure  8.  Simulated  SWIFT/486  results  and  analytical  solutions  for  a  radial 
solute  transport  problem  with  dispersivity  coefficient  of  0.15  m 


Problem  5,  Transport  from  a  continuous  point  source  in  a  uniform 
two-dimensional  flow  field 

Objectives.  The  purpose  of  Example  Problem  5  is  to  evaluate  SWIFT/486 
for  transport  options  solving  a  continuous  source  in  a  Cartesian  coordinate 
system  and  to  check  decay  and  adsorption/desorption  options. 
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Figure  9.  Simulated  SWIFT  results  and  analytical  solutions  for  a  radial 

solute  transport  problem  with  dispersivity  coefficient  of  0.015  m 


Problem  statement.  In  this  problem,  a  conservative  tracer  is  injected  into 
an  aquifer  via  a  fully  penetrating  well  (Beljin  1991,  1993).  A  plume  develops 
from  the  source  and  spreads  out  to  the  sides  of  the  aquifer.  If  one  assumes 
uniform  mixing  of  mass  throughout  the  aquifer  thickness,  two-dimensional 
(x-y  plane)  tracer  solute  transport  can  be  described  as  follows: 


D  ^ 


V~  -  IRC  ^  ^  =  R  — 


dx 


4> 


dt 


(10) 


C(^y,0)  =  0 
Qc^yJ)  =  QC^SQcy) 
C  (±  00,  ±  00,  r )  =  0 


(11) 


where  X  is  a  first-order  decay  constant  [1/T];  is  the  retardation  coefficient 
[dimensionless];  is  the  mass  injection  rate  of  solute  per  unit  volume  of 
aquifer  [M  T‘  L‘^];  Q  is  the  volumetric  injection  rate  of  fluid  per  unit  of  aqui¬ 
fer  thickness  [L^  T  'j;  Q  is  the  concentration  of  the  injected  fluid  [M  L‘^]; 
and  dix.y)  is  the  Dirac  delta  function  [1/L^].  In  SWIFT/486,  the  retardation 
coefficient  is  calculated  using  the  formation  density  (BWRN),  the  adsorp¬ 
tion  coefficient  Kj  (DIS),  and  the  formation  porosity  <t>  (POROS)  as: 

«=  1  ♦  (12) 


Chapter  3  Model  Performance 


21 


Note  that  BWRN,  DIS,  POROS  are  variable  names  used  in  SWIFT/486  for 
formation  density,  adsorption  coefficient,  and  porosity,  respectively.  For 
retardation  factor  of  7?  =  2,  the  adsorption  coefficient  was  calculated  to  be: 

>^</  =  1  +  -^3^^  X  1690  =  3.263  x  lO  '*  (13) 


SWIFT/486  calculates  the  concentration  in  mass  fraction  (normalized).  To 
convert  mass  fraction  to  other  concentration  units,  for  example  to  milligrams 
per  liter,  mg/f,  the  following  relationship  can  be  used: 

^  =  mass  fraction  ^  fluid  density  (kglm^)  x  i,000  (14) 


where  1 ,000  is  a  conversion  factor. 

Input.  The  input  parameters  used  in  this  example  are  from  a  measurement 
at  a  field  site  located  in  Long  Island,  New  York  (Finder  1973,  Wilson  and 
Miller  1978,  Beljin  1988).  The  original  data  were  measured  for  a  plume  of 
hexavalent  chromium.  However,  for  this  example,  it  was  assumed  that  the 
chemical  did  not  change  chemically;  therefore,  the  type  of  chemical  is  imma¬ 
terial.  The  major  input  data  are  given  in  Table  9.  Note  that  flow  velocity 
was  input  directly  into  the  analytical  solution  while  for  the  numerical  solution 
it  was  calculated  using  the  hydraulic  conductivity  and  gradient. 

Numerical  specification.  A  Cartesian  coordinate  system  with  32  elements 
in  the  x-direction,  5  elements  in  the  y-direction,  and  1  element  in  the  z- 
direction  was  used.  Spacing  in  the  x-direction  (DX)  was  set  to  60  m,  in  the 
y-direction  (DY)  was  set  to  30  m,  50  m,  and  60  m,  respectively,  for  the  three 
different  simulations,  and  in  the  z-direction  (DZ)  was  set  to  33.5  m  (depth  of 
aquifer).  The  radioactive  transport  option  of  SWIFT/486  was  used  to  simulate 
a  nonradioactive  chemical  with  and  without  adsorption  in  effect. 

Output.  The  output  for  Example  Problem  5  is  included  in  Appendix  A. 
Concentration  versus  the  distance  profiles  are  given  in  Figures  10-15. 

Results.  The  results  of  six  different  scenarios  are  shown  in  Figures  10-15. 
In  Figures  10,  11,  and  12,  all  the  input  parameters  including  Ax  and  Az  were 
fixed,  but  Ay  was  changed  to  30,  50,  and  60  m.  As  shown  in  Figures  10,  11, 
and  12,  the  variation  of  the  grid  spacing  in  the  y-direction  had  little  effect  on ' 
the  results.  In  another  scenario,  the  retardation  coefficient  was  changed  from 
R  =  1  to  7?  =  2  (Figures  12  and  15),  and  as  expected,  the  solute  transport 
delayed.  For  retardation  coefficient  of  7?  =  1  (no  adsorption),  Ar  =  60  m. 

Ay  =  60  m,  and  A/  =  100  days,  the  simulation  results  at  the  times  of 
1,000  days,  2,000  days,  and  2,800  days  (Figures  12,  13,  and  14)  showed  a 
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Table  9 

Input  Parameter  for  Example  Problem  5  (Beljin  Cartesian  Problem) 

Parameter 

Value 

Darcy  velocity,  v 

0.161  m/d  (0.525  ft/d) 

Seepage  velocity,  U 

0.460  m/d  (1.500  ft/d) 

Porosity,  n 

0.35 

Longitudinal  dispersivity, 

21.3  m  (69.9  ft) 

Transverse  dispersivity,  Oj 

4.27  m  (14.0  ft) 

Aquifer  saturated  thickness,  b 

33.5  m  (110  ft) 

Point  source  strength,  Q 

23.59  kg/d  (52  Ib/d) 

oq> 

704.0  g/(m  d) 

Time,  t 

1000,  2000,  and  2800  d 

Case  1 ,  retardation  factor,  R 

1.0 

Decay  constant  ,  X 

0.0  1/d 

Case  2,  retardation  factor,  R 

2,0 

Decay  constant,  X 

0.00019  1/d 

40 

g  ^0 

c 

o 

O  20 
1  0 
0 


DX  =  60m 
DY  =  30m 
DT  =  100  days 
Time  =  1000  days 
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Figure  10,  Simulated  SWIFT  results  and  analytical  solutions  for  two- 

dimensional  continuous  source  problem  with  DY  =  30  m  and 
retardation  coefficient  /?  =  1 
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Figure  1 1 .  Simulated  SWIFT  results  and  analytical  solutions  for  two- 

dimensional  continuous  source  problem  with  DY  =  50  m  and 
retardation  coefficient  R  =  1 


Distance  (m) 


Figure  12.  Simulated  SWIFT  results  and  analytical  solutions  for  two- 

dimensional  continuous  source  problem  with  DY  =  60  m  and 
retardation  coefficient  R  =  1 


24 


Chapter  3  Model  Performance 


Analytic  Qi 


DX  =  60  m 
DY  =  60  m 
DT  =  1 00  days 
Time  =  2000  days 


200  400  600  800  1000  1200  1400  1600  1800  2000 


Distance  (m) 


Figure  13.  Simulated  SWIFT  results  and  analytical  solutions  for  two- 

dimensional  continuous  source  problem  at  time  =  2,000  days 
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Figure  14.  Simulated  SWIFT  results  and  analytical  solutions  for  two- 

dimensional  continuous  source  problem  at  time  =  2,800  days 
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Figure  15.  Simulated  SWIFT  results  and  analytical  solutions  for  two- 

dimensional  continuous  point  source  problem  with  DY  =  60  and 
retardation  coefficient,  R  =  2 


good  agreement  with  the  analytical  solution.  Beljin  (1988)  compared  the  same 
analytical  solution  with  models  such  as  MOC,  SEFTRN,  and  RANDOM 
WALK  and  reported  a  similar  conclusion. 


Problem  6,  Transport  of  a  solute  slug  in  a  uniform  groundwater  flow 
field 

Objectives.  The  purpose  of  this  test  is  to  evaluate  SWIFT/486  for  slug 
transport  of  a  conservative  solute  using  the  radioactive  transport  option  of  the 
code  in  a  Cartesian  coordinate  system. 

Problem  statement.  This  problem  considers  the  slug  injection  of  a 
conservative  solute  into  a  uniform  two-dimensional  flow  field.  The  difference 
between  this  example  and  Example  Problem  5  is  the  initial  source  condition. 
The  initial  concentration  in  this  problem  is  a  function  of  the  source  mass  (M) 
as: 


C  {x,y,Qi)  =  —  8ix,y) 
n 


where 

M  (M/L)  =  mass  of  the  solute  injected  instantaneously  into  the 
aquifer  per  unit  length  of  aquifer  thickness 

n  =  porosity 
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5(x,y)  (1/L"]  =  Dirac  delta  function 

Input.  Input  parameters  for  Example  Problem  6  (Table  10)  are  from 
Beljin  (1988).  SWIFT/486  requires  input  concentrations  only  in  units  of  mass 
fraction.  For  this  example,  the  aquifer  thickness  is  33.5  m,  and  the  total  mass 
injected  over  the  aquifer  thickness  is  117.25  kg  (3.5  kg/m  x  33.5  m).  The 
mass  of  pore  water  for  a  cell  is  Ax  Ay  Az  </•  where  Ax,  Ay,  and  Az  are  the 
grid  spacing  in  the  x-,  y-,  and  z-directions,  respectively,  </>  is  the  porosity,  and 
is  the  water  density.  The  mass  fraction  then  is  defined  as  follows; 


r  ^  —  total  mass  injected  /\£\ 

mass  fraction  =  - ^ ; -  (to/ 

mass  of  pore  water  m  column 


Table  10 

Input  Parameters  for  Example  Problem  6 

Parameter 

Value 

Darcy  velocity,  v 

2.0  m/d  (6.56  ft/d) 

Seepage  velocity,  v 

5.71  m/d  (18.75  ft/d) 

Porosity,  n 

0.35 

Longitudinal  dispersivity,a^ 

4.0  m  (13.12  ft) 

Transverse  dispersivity,a;^ 

1.0  m  (3.28  ft) 

Retardation  factor,  R 

1 

Decay  constant,  A. 

0.0  1  /d 

Solute  mass  per  unit  aquifer 
thickness 

3.5  kg/m  (2.35  Ib/ft) 

! 

Mass  fraction 

0.0004 

Time,  t 

3.96,  10.59,  and  16.59  d 

Numerical  specification.  Two  different  rectangular  grids,  one  with 
40  elements  in  the  x-direction,  5  elements  in  the  y-direction,  and  1  element  in 
the  z-direction,  and  another  with  the  same  number  of  elements  in  the  x-  and  z- 
directions  but  19  elements  in  the  y-direction  were  used.  The  radionuclide 
transport  option  was  used  to  solve  this  problem.  The  seepage  velocity  or 
Darcy  velocity  required  by  the  analytical  solution  cannot  be  entered  directly 
into  SWIFT/486.  The  velocity  is  input  implicitly  into  the  numerical  solution 
by  assigning  estimated  pressure  heads  at  the  boundaries  and  the  calculated 
hydraulic  conductivity  of  the  aquifer.  Other  specifications  are  given  in  the 
figures. 
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Output.  Output  for  Example  Problem  6  is  included  in  Appendix  A.  Con¬ 
centration  (parts  per  million)  versus  distance  (meter)  profiles  are  plotted  in 
Figures  16-21. 

Results.  The  results  of  six  different  scenarios  are  given  in  Figures  16,  17, 
18,  20,  21,  and  22.  In  Figures  16,  17,  and  18,  although  the  locations  of  peak 
concentration  for  the  numerical  results  and  analytical  solution  are  the  same, 
there  are  some  differences  between  the  numerical  and  analytical  results  around 
the  peak  concentrations.  By  increasing  the  domain  size  in  the  y-direction 
from  25  m  to  95  m,  better  results  were  obtained  as  shown  in  Figures  22  and 
23.  The  reason  for  obtaining  better  results  for  domain  with  y  =  96  m  is  due 
to  providing  enough  space  in  y-direction  for  plume  to  spread  out  completely. 
The  lateral  increase  of  domain  size  did  not  change  the  results  for  time  =  3.96 
days.  The  discrepancy  between  the  numerical  and  analytical  results  at  time 
3.96  days  is  because  the  software  package  of  analytical  solutions  does  not 
produce  accurate  results  at  this  time  (  Beljin  1988).  Beljin  (1988)  compared 
MOC,  SEFTRAN,  and  RANDOM  WALK  simulated  results  with  the  same 
analytical  solution  and  obtained  similar  agreements  between  the  solutions. 

SWIFT/486  is  a  relatively  fast  code.  To  provide  the  reader  with  an  idea 
on  the  running  speed  of  the  code,  the  estimated  CPU  times  for  this  problem 
are  given  in  Figures  19  and  23.  Note  that  the  PC  machine  used  in  these  cal¬ 
culations  was  a  486/DX2,  66  MHZ  with  24  megabyte  RAM.  As  expected,  a 
comparison  between  Figures  19  and  23  shows  that  more  CPU  time  is  used  for 
a  case  with  19  elements  in  the  y-direction  than  for  5  elements  in  the 
y-direction. 

This  test  showed  that  the  radioactive  transport  option  of  SWIFT/486  can  be 
used  to  simulate  a  conservative  solute  transport  correctly  within  the  assump¬ 
tions  of  the  selected  analytical  solution. 
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Figure  16.  Simulated  SWIFT/486  results  and  analytical  solution  for  two 
dimensional  slug  transport  of  a  solute  at  time  =  3.96  days 
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Figure  17.  Simulated  SWIFT/486  results  and  analytical  solution  for  two 
dimensional  slug  transport  of  a  solute  at  time  =  10.59  days 
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Figure  18.  Simulated  SWIFT/486  results  and  analytical  solution  for  two- 
dimensional  slug  transport  of  a  solute  at  time  =  16.59  days 


Figure  1 9.  SWIFT/486  CPU  time  calculated  for  the  simulations  with 
NX  =  40  and  NY  =  5 
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Figure  20.  Simulated  SWIFT/486  results  and  analytical  solution  for  two 
dimensional  slug  transport  of  a  solute  with  NX  =  40  and 
NY  =  19  at  time  =  3.96  days 
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Figure  21 .  Simulated  SWIFT/486  results  and  analytical  solution  for  two 
dimensional  slug  transport  of  a  solute  with  NX  =  40  and 
NY  =  19  at  time  =  10.59  days 


Chapter  3  Model  Performance 


20 


OX  =  5  m 
DY  =  5  IT) 

DT  =  1  hour 
Time  =  1  6.59  days 
NX  =  40 
NY  =  1  9 


20  40  60  80  100  120  140  160  180  200 

Distance  in  X'Direction  (m) 


Figure  22.  Simulated  SWIFT/486  results  and  analytical  solution  for  two- 
dimensional  slug  transport  of  a  solute  with  NX  =  40  and 
NY  =  19  at  time  =  16.59  days 


Figure  23.  SWIFT/486  CPU  time  calculated  for  simulations  with  NX 
and  NY  =  19 
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4  Conclusions 


Different  versions  of  SWIFT/486  have  been  evaluated,  tested,  and  applied 
at  several  sites  by  petroleum  engineers  and  hydrogeologists  (Ward  et  al. 

1987).  A  list  of  some  of  the  SWIFT/486  applications  is  given  in  Table  1 1 
(GeoTrans  1994).  The  evaluation  performed  here  complements  previous 
SWIFT  evaluations  and  applications  (Table  11).  In  this  investigation,  the 
numerical  accuracy  of  SWIFT/486  has  been  evaluated  by  a  comparison  of 
simulation  results  with  analytical  solutions.  The  selected  analytical  problems 
include  a  variety  of  initial  and  boundary  conditions.  The  model  also  was 
reviewed  for  efficiency  of  coding  by  checking  its  speed  of  run,  convenience  of 
input/output  by  checking  data  input  and  output  information,  program  porta¬ 
bility  by  running  it  on  two  different  computer  systems  (DOS  and  UNIX),  and 
available  diagnostic  messages  received  during  the  simulations.  Overall, 
SWIFT/486  is  a  relatively  efficient  code,  requires  an  optimal  amount  of  com¬ 
puter  storage,  and  has  sufficient  diagnostic  flags.  SWIFT/486  simulations 
closely  matched  the  analytical  solutions  to  several  simplified  problems. 
SWIFT/486,  however,  has  a  few  deficiencies  which  make  its  use  inconvenient 
as  described  below. 

The  lack  of  a  preprocessor  for  preparing  input  files  is  one  of  the  draw¬ 
backs  which  impedes  many  users  from  selecting  SWIFT/486.  SWIFT/486 
requires  extensive  data  and  input  parameters.  Input  files  require  formatted 
structure.  A  SWIFT/486  user  should  also  have  adequate  knowledge  of  ther¬ 
modynamics,  rock  mechanics,  and  groundwater  hydrology.  There  are  many 
options  in  SWIFT/486  which  make  the  application  of  the  code  very  flexible 
for  an  experienced  user  and  difficult  for  novices.  Some  of  the  options  are 
described  in  this  report;  these  and  others  can  be  found  in  Ward  et  al.  (1993). 
One  of  the  useful  options  in  SWIFT/486,  for  instance,  is  that  the  daily  varia¬ 
tions  in  the  injected  density  can  be  assigned  by  varying  the  solute  mass  frac¬ 
tion.  There  are  some  options  in  which  the  user  must  manipulate  the  input 
data  in  order  to  run  the  code.  For  example,  SWIFT/486  cannot  be  used  for 
pure  advection  problems  (e.g.,  dispersivity  coefficient  =  0).  For  problems 
with  zero  dispersivity,  values  near  zero  must  be  input;  otherwise,  the  control 
statement  in  the  code  will  stop  the  simulation.  This  minor  change  of  input 
parameter  will  allow  the  user  to  run  the  code  and  at  the  same  time  not  deviate 
much  from  the  actual  input  parameter.  There  are  some  options  available  in 
SWIFT/486  which  are  not  defined  in  the  user’s  manual.  For  example,  the 


Chapter  4  Conclusions 


33 


Table  1 1 

SWIFT  Model  Applications  (GeoTrans  1994) 

Site 

State 

Dimensions 

Flow 

Transport 

Hazardous 

Radioactive 

Salt  Water 

Chem-Dyne 

OH 

3 

Y 

Y 

Ottatl  &  Goss 

NH 

3 

Y 

Y 

Woburn 

MA 

3 

Y 

Y 

Fernald 

OH 

3 

Y 

Y 

Power  Rd 

OH 

3 

Y 

Y 

S-Area 

NY 

3 

Y 

Y 

Savannah 

SC 

3 

Y 

Y 

Conserv  Chem 

MO 

3 

Y 

Y 

Confidential- 1 

CA 

3 

Y 

Y 

Confidential-2 

CN 

2 

Y 

Y 

Confidentia!-3 

FL 

3 

Y 

Y 

Y 

Confidential-4 

Ml 

3 

Y 

Y 

Confidential-5 

OK 

3 

Y 

Y 

SWF  WAT 

MGMN  DIS 

FL 

2 

Y 

BWIP 

WA 

3 

Y 

Y 

WIPP 

NM 

3 

Y 

Y 

Volusia,  County 

FL 

3 

Y 

Y 

Saudi  Arabia 

3 

Y 

Y 

Borden  Landfill 

CAN 

3 

Y 

Y 

Babylon  Landfill 

NY 

2 

Y 

Y 

Ates  Mobile 

AL 

3 

Y 

Musquodoboit 

CAN 

2 

Y 

East  Kent  Chalk 

UK 

2 

Y 

Sacramento 

CA 

3 

Y 

Y 

Confidential 

AK 

3 

Y 

Y 

parameter  METHOD  in  M3-1  record  of  input  file  can  have  a  zero  value;  how¬ 
ever,  this  is  not  defined  in  the  user’s  manual.  The  future  user’s  manual 
should  cover  all  the  new  options  available  in  this  version  of  SWIFT/486. 
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SWIFT/486  has  been  developed  to  solve  saturated  zone  problems;  how¬ 
ever,  there  is  an  option  in  SWIFT/486  that  allows  the  user  to  simulate  the 
dewatering  (draining)  of  a  saturated  zone  where  the  saturation  may  vary  from 
1.0  to  0.001.  This  option  in  SWIFT/486  is  described  under  the  Free  Water 
Surface  option  (Reeves  et  al.  1986a).  As  described  by  Reeves  et  al.  (1986a), 
there  are  some  problems  with  this  option.  Therefore,  this  feature  of 
SWIFT/486  was  not  evaluated  here.  SWIFT/486  is  categorized  as  a  saturated 
zone  model. 

SWIFT/486  has  been  developed  using  the  parameters  and  terms  derived  for 
petroleum  engineers.  For  example,  it  calculates  pressure  in  PSI  units  rather 
than  hydraulic  head.  Furthermore,  the  aquifer  storage  coefficient  and  the 
aquitard  specific  coefficient  are  calculated  from  the  input  values  for  water 
compressibility,  rock  compressibility,  and  effective  porosity.  For  hydro¬ 
geology  projects,  the  unit  of  the  calculated  parameters  by  SWIFT/486  could 
be  converted  into  the  units  that  are  used  by  hydrogeologists,  civil  engineers, 
or  environmental  engineers. 

For  a  radial  (cylindrical  geometry)  coordinate  system,  SWIFT/486  calcu¬ 
lates  the  Peclet  number  wrongly,  therefore,  the  warning  received  by  the  user 
is  not  based  on  actual  Peclet  number  of  the  problem  and  should  be  ignored.’ 
The  user  may  wish  to  calculate  the  Peclet  number  for  this  type  of  problem 
using  other  available  formulation  (Bear  1979). 

Overall,  this  investigation  showed  that  for  selected  problems  with  simplify¬ 
ing  assumptions,  SWIFT/486  performed  very  well. 


*  Personal  Communication,  1994.  D.  Ward,  Vice  President,  GeoTrans,  Inc.,  Sterling,  VA  20166. 
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Appendix  A 

Input  and  Output 

of  SWIFT/486  Simulations 


Appendix  A  provides  a  complete  list  of  input  and  output  parameters 
used  in  this  evaluation. 
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INDEX  OF  RESERVOIR  HETEROGENEITY  ....  HTC 

NO  OF  RADIOACTIVE  COMPONENTS . NCP  .. 

NUMBER  OF  ROCK  TYPES . NRT  ..  1 
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TEMPERATURE-VISCOSITY  TABLE 
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CYLINDRICAL  GLOBAL  SYSTEM  DATA  *** 

WELLBORE  RADIUS . RW  ....  0.1143  (M) 

RADIUS  TO  CENTER  OF  FIRST  COLUMN  ...  R1  ....  0.4755  (M) 

RESERVOIR  EXTERIOR  RADIUS . RE  ....  6096.  (M) 

DEPTH  TO  CENTROID  OF  BLOCK  (1,1,1) .  DEPTH  .  O.OOOOE+00 
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DATA  FOR  CARTER-TRACY  WATER  INFLUX  CALCULATIONS 
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DEPTH  OF  BLOCK  CENTERS  BELOW  REFERENCE  PLANE  (M) 
(Measured  positive  downwards) 
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♦**  STATE  VARIABLE  INITIALIZATION 


♦♦♦  RECURRENT  DATA  SPECIHCATION  BEGINNING  AT  TIME  =  O.OOOOE+00  (SECS) 
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TIME  STEPPING  AND  OUTPUT  CONTROL  OPTIONS 
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AVERAGE  PRESSURE  5.9770E+04  (PA)  HEAT  LOSS  TO  OVER/UNDRBRDN  O.OOOOE+00  (J) 


WELL  OPERATION  SUMMARY 
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34  5.9770E+04  21.100  0.0000  35  5.9770E-H04  21.100  0.0000  36  5.9770E-h04  21.100  0.0000 

37  5.9770E+04  21.100  0.0000  38  5.9770E-I-04  21.100  0.0000  39  5.9770E+04  21.100  0.0000 

40  5.9770E+04  21.100  0.0000  41  5.9770E+04  21.100  0.0000  42  5.9770E+04  21.100  0.0000 
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GLOBAL  (FRACTURE)  DEPENDENT  VALUES 
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TIME  STEP  NUMBER  3  NUMBER  OF  OUTER  ITERATIONS  1  CURRENT  TIME  STEP  1.2960E+05  SECS 
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ELAPSED  SIMULATION  TIME  7.0286E+05  SECS  (  8.135  DAYS  ,  2.2288E-02  YEARS  ) 
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43  5,8443E+04  21.100  0,0000  44  5.8904E+04  21.100  0.0000  45  5.9255E+04  21.100  0.0000 

46  5.9497E+04  21.100  0.0000  47  5,9645E+04  21.100  0.0000  48  5.9722E+04  21.100  0  0000 

49  5.9755E+04  21.100  0.0000  50  5.9767E+04  21.100  0,0000 
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28  4.5119E+04  21.100  0.0000  29  4.6012E+04  21.100  0.0000  30  4,6904E+04  21.100  0,0000 

31  4.7796E+04  21.100  0.0000  32  4.8686E+04  21.100  0.0000  33  4,9574E+04  21.100  0.0000 
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22  3.8718E+04  21.100  0.0000  23  3.9612E+04  21.100  0.0000  24  4.0506E+04  21.100  0.0000 
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10  2.6983E+04  21.100  0.0000  11  2.7877E+04  21.100  0.0000  12  2.8771E+04  21.100  0.0000 

13  2.9665E+04  21.100  0.0000  14  3.0559E+04  21.100  0.0000  15  3.1453E+04  21.100  0.0000 
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OBSERVATION  WELL  (R=65.61  FT  OR  20.0  M) 
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PROB.  2.0  (MR)  ++  FLOW  VERIRCATION  -  SI  (METRIC)  -  RADIAL  COORDS 
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TEMPERATURE-VISCOSITY  TABLE 
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DATA  FOR  CARTER-TRACY  WATER  INFLUX  CALCULATIONS 
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AQUIFER  COEFFICIENTS  (DIMENSIONLESS) 


Appendix  A  SWIFT/486  Simulations 


AQUIFER-INFLUENCE  FUNCTION  BLOCK  NUMBERS 
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DEPTH  OF  BLOCK  CENTERS  BELOW  REFERENCE  PLANE  (M) 
(Measured  positive  downwards) 
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GLOBAL  PORE  VOLUME  (M**3) 
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GLOBAL  ROCK  TYPES 
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INITIAL  GLOBAL  PRESSURE  AT  ELEVATION  H  (PA) 
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*  STATE  VARIABLE  INITIALIZATION  **♦ 
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(GLOBAL+LOCAL)  MASS  OR  HEAT  BALANCE  1.0000  *♦*♦♦** 

WELL  SUMMARY 

TOTAL  PRODUCTION  121.0  1.0695E+07  O.OOOOE+00 

TOTAL  INJECTION  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
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GLOBAL  PRESSURE  AT  ELEVATION  (PA) 


GLOBAL  TEMPERATURE  (DEG.C) 
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GLOBAL  (FRACTURE)  DEPENDENT  VALUES 
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GLOBAL  PRESSURE  AT  ELEVATION  (PA) 
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GLOBAL  PRESSURE  AT  DATUM  (PA) 
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GLOBAL  TEMPERATURE  (DEG.C) 
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GLOBAL  (FRACTURE)  DEPENDENT  VALUES 
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GLOBAL  PRESSURE  AT  ELEVATION  (PA) 
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GLOBAL  PRESSURE  AT  DATUM  (PA) 
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GLOBAL  TEMPERATURE  (DEG.C) 
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AQUIFER  INFLUX  RATES  (POSlTIVE-IN  :  NEGATIVE-OUT) 
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X-DIRECnON  DISTANCE  TO  GRID  BLOCK  CENTER  (M) 
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**♦  STATE  VARIABLE  INITIALIZATION 
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»♦*  RECURRENT  DATA  SPECIHCATION  BEGINNING  AT  TIME  =  O.OOOOE+00  (SECS) 
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GLOBAL  TEMPERATURE  (DEG.C) 
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INDEX  OF  RESERVOIR  HETEROGENEITY  ....  HTC 

NO  OF  RADIOACTIVE  COMPONENTS . NCP  .. 

NUMBER  OF  ROCK  TYPES . NRT ..  1 
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TEMPERATURE-VISCOSITY  TABLE 
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***  CYLINDRICAL  GLOBAL  SYSTEM  DATA  ♦♦♦ 

WELLBORE  RADIUS . RW....  0.1 143  (M) 

RADIUS  TO  CENTER  OF  FIRST  COLUMN  ...  R1  ....  0.4755  (M) 

RESERVOIR  EXTERIOR  RADIUS . RE....  50.00  (M) 

DEPTH  TO  CENTROID  OF  BLOCK  (1,1,1) .  DEPTH  .  O.OOOOE+00 
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DEPTH  OF  BLOCK  CENTERS  BELOW  REFERENCE  PLANE  (M) 
(Measured  positive  downwards) 
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AQUIFER  INFLUX  RATES  (POSITIVE-IN  :  NEGATIVE-OUT) 
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SANDIA  Waste-rsolation 
Flow  and  Transport  in 
Porous  and/or  Fractured  Media 
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Copyright  GeoTrans,  Inc.  1993 
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NUMBER  OF  BLOCKS  IN  X-DIRECTION . NX  . 

NUMBER  OF  BLOCKS  IN  Y-DIRECTION . NY  ., 

NUMBER  OF  BLOCKS  IN  Z-DIRECTION . NZ  .. 
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GLOBAL  SYSTEM  GRIDDING  **♦ 
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X-DIRECnON  DISTANCE  TO  GRID  BLOCK  CENTER  (M) 
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♦**  SPECinCATION  OF  HOMOGENEOUS  GLOBAL  SYSTEM 
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DEPTH  OF  BLOCK  CENTERS  BELOW  REFERENCE  PLANE  (M) 
(Measured  positive  downwards) 


GLOBAL  PORE  VOLUME  (M**3) 


Appendix  A  SWIFT/486  Simulations 


A207 


All  values  for  this  array  equal  6.3650E-05 
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***  RECURRENT  DATA  SPECIFICATION  BEGINNING  AT  TIME  =  1.7280E+07  (SECS) 
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NORMAL  TERMINATION  (mME=  28  ;  TIME  =  2.4192E+08 ) 
CPU  elapsed  time  =  10.11  seconds 
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SANDIA  Waste-Isolation 
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INDEX  OF  RESERVOIR  HETEROGENEITY  ....  HTG 

NO  OF  RADIOACTIVE  COMPONENTS . NCP .. 

NUMBER  OF  ROCK  TYPES . NRT ..  1 


00 

S  o 
o  ■ 

1^4  H  ® 
•  S  HJ 
:  S  U 

:  :  OC 


Z  o 
tq  < 
J  2  O 
S  J  Q 

^  9  < 

DC 


fc  5  X 
5e2  < 
o  £ 


:  O  X 

'  ^  S 
X  a: 

fflZo 

<. 

^  CO  Q 
i  X  O 

gSs 

y  m  s 

2g  i 

f?  <  ^ 
^  X  i 

iij 

i§g 

£=:  tu  H 
^ 

CXd  o 

<  w  ^ 
u-  }ir 

O  O  o 
O  O  S 
z.  z.  Q 

ssl 

S  S  US 


<  H  U  9 

S  S  9  K 
S  2  ft 
<zo8 

^  O  a:  lZ 
Q  LG  X 


A228 


Appendix  A  SWIFT/486  Simulations 


NO  MASS  COMPONENT  PARENT  COMP  DECAY  ADJUST  FREUNDLICH  COEFFS.  KAPPA  (CU.M/KG)  ♦♦ETA  AND  ETA 
NO  FRACTION  HLIFE(YRS)  LAMDA  ROCK  TYPES 
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TEMPERATURE-VISCOSITY  TABLE 
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GLOBAL  SYSTEM  GRIDDING 
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X-DIRECTION  DISTANCE  TO  GRID  BLOCK  CENTER  (M) 
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Y-DIRECTION  DISTANCE  TO  LEADING  BLOCK  EDGE  (M) 
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SPECinCATION  OF  HOMOGENEOUS  GLOBAL  SYSTEM 
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DEPTH  OF  BLOCK  CENTERS  BELOW  REFERENCE  PLANE  (M) 
(Measured  positive  downwards) 
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***  SALT  DISSOLUTION  *** 
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GLOBAL  PORE  VOLUME  (M**3) 
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All  values  for  this  array  equal  7.7050E-04 
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INITIAL  GLOBAL  PRESSURE  AT  DATUM  ELEVATION  (PA) 
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o  ^ 
o 

o 


C'l  , 
•  , 
»r»  00 


O 
Z 

CQ 

W  _ 


w  S 

-:pg 


tu  p  5  y 

J  J 

K  CQ  £  ^ 


Appendix  A  SWIFT/486  Simulations 


♦  RECURRENT  DATA  SPECIFICATION  BEGINNING  AT  TIME  =  7200.  (SECS) 


A262 


Appendix  A  SWIFT/486  Simulations 


GLOBAL  PRESSURE  AT  ELEVATION  (PA) 


o 


On 


OO 


t-* 


VO 


Tt 


m 


VD 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

VD 

VO 

o 

o 

O 

O 

o 

o 

O 

o 

O 

o 

+ 

+ 

4" 

4- 

4* 

4- 

+ 

4- 

+ 

4* 

(2 

u 

U 

u 

m 

PCI 

PQ 

PQ 

PJ 

03 

r- 

<s 

cs 

CM 

CM 

CM 

«n 

m 

r4 

cs 

CM 

CM 

CM 

«n 

in 

m 

in 

in 

in 

m 

On 

On 

Tj- 

ON 

On 

Tf 

ON 

Tf 

s 

8 

8 

8 

o 

p 

OH 

1-^ 

HH 

hH 

HH 

HH 

HH 

hh 

NO 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

o 

o 

O 

o 

o 

O 

o 

o 

O 

o 

+ 

+ 

4- 

4- 

4" 

4* 

4- 

4* 

4* 

4- 

UJ 

(D 

in 

U 

Uj 

PJ 

PJ 

PJ 

03 

o 

o 

O 

o 

o 

m 

m 

in 

in 

m 

VO 

VO 

VO 

VO 

VO 

r4 

cs 

CM 

CM 

CM 

5 

5 

5 

5 

5 

Tf 

in 

Tf 

in 

Tf 

m 

Tf 

in 

Tf 

m 

•n 

m 

m 

in 

in 

p 

p 

p 

p 

p 

OH 

HH 

HH 

*.H 

HH 

VO 

VO 

NO 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

O 

o 

o 

o 

o 

O 

O 

o 

o 

o 

+ 

4" 

4" 

4- 

4- 

4- 

4* 

4- 

4- 

4- 

m 

UJ 

U 

UJ 

U 

PO 

PQ 

PQ 

U 

03 

Tf 

rj- 

Tf 

Tf 

Tf 

On 

ON 

C3N 

ON 

ON 

VO 

VD 

VO 

VO 

VO 

CnJ 

<s 

CM 

CM 

CM 

m 

m 

m 

cn 

CO 

CO 

CO 

CO 

CO 

CO 

ON 

On 

On 

ON 

ON 

O 

o 

o 

O 

o 

V-) 

»n 

in 

*n 

J 

OH 

OH 

hH 

*-H 

HH 

hH 

HH 

HH 

HH 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

VD 

VO 

o 

o 

o 

O 

o 

o 

o 

O 

O 

o 

+ 

+ 

4" 

4- 

4* 

4- 

4- 

4- 

4" 

4- 

U 

ua 

PJ 

U 

PJ 

PCI 

va 

03 

t-* 

r- 

r- 

C4 

cs 

CM 

CM 

CM 

VO 

VO 

VO 

VO 

VO 

CO 

CO 

CO 

CO 

CO 

CN 

<N 

cs 

fS 

CM 

<N 

CM 

CM 

CM 

CM 

Tf- 

Tf 

Tf 

Tf 

m 

in 

m 

in 

in 

VO 

vq 

vq 

vq 

vq 

hh 

oJ 

OH 

H-i 

hH 

HH 

HH 

HH 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

o 

O 

O 

o 

o 

cs 

O 

o 

O 

o 

O 

+ 

+ 

4- 

+ 

4- 

4- 

4- 

4* 

4- 

4- 

U 

Ui 

tu 

u 

m 

PJ 

PLI 

03 

VO 

VO 

VO 

VO 

VO 

r- 

t> 

r- 

r- 

CO 

CO 

CO 

CO 

CO 

On 

On 

ON 

On 

O'. 

On 

o 

o 

o 

O 

o 

VO 

VO 

vq 

vq 

vq 

r4 

CM 

CM 

CM 

CM 

OH 

o-^ 

hh 

HH 

HH 

hH 

hh 

HH 

hh 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

O 

O 

o 

O 

o 

O 

O 

o 

o 

O 

+ 

+ 

4" 

4- 

+ 

4* 

4- 

4- 

4* 

4- 

u 

IXI 

m 

00 

PQ 

PJ 

U 

PJ 

03 

Tj- 

Tt 

Tf 

*vt 

Tf 

ON 

ON 

On 

On 

On 

CO 

CO 

CO 

CO 

CO 

o 

o 

o 

o 

o 

o 

o 

O 

O 

O 

Tt 

Tf 

Tf 

Tf 

in 

in 

in 

m 

in 

t-- 

cs 

cs 

CM 

CM 

CM 

OH 

o-^ 

HH 

HH 

HH 

HH 

HH 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

o 

O 

o 

O 

o 

O 

o 

O 

O 

O 

+ 

+ 

4- 

+ 

4- 

4- 

4* 

4- 

4- 

4* 

U 

ui 

tl] 

PO 

ta 

p:^ 

PJ 

03 

OO 

00 

OO 

00 

OO 

CO 

CO 

CO 

CO 

CO 

r- 

r- 

r-- 

VO 

Tf 

Tf 

Tf 

Tf 

On 

ON 

On 

On 

On 

On 

On 

On 

On 

On 

OO 

OO 

00 

OO 

OO 

ON 

On 

ON 

On 

On 

r- 

r- 

r-; 

CNl 

CM 

CM 

CM 

CM 

oH 

o-.^^ 

HH 

HH 

HH 

hh 

HH 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

O 

o 

O 

o 

o 

o 

O 

O 

o 

O 

+ 

+ 

4* 

+ 

4- 

4- 

4* 

4- 

4- 

4- 

U 

u 

W 

(D 

U 

\Xi 

ta 

m 

03 

VO 

VO 

VO 

VO 

VO 

OO 

00 

00 

00 

OO 

Tf 

Tf 

Tf 

Tf 

Tf 

OO 

OO 

OO 

00 

00 

OO 

OO 

OO 

OO 

OO 

m 

CO 

CO 

CO 

Tf 

Tf 

Tf 

Tf 

Tf 

OO 

00 

00 

OO 

OO 

Tf 

CO 

p 

p 

p 

p 

OH 

OO^ 

HH 

HH 

HH 

VO 

VO 

VO 

NO 

VO 

VO 

VO 

VO 

NO 

VO 

O 

o 

O 

o 

o 

O 

o 

O 

o 

o 

+ 

+ 

+ 

4" 

4- 

4* 

4- 

4- 

4" 

4- 

U 

UJ 

tu 

va 

CO 

to 

PJ 

PJ 

PJ 

03 

m 

«r> 

in 

m 

O 

O 

O 

O 

o 

OO 

OO 

00 

OO 

00 

in 

m 

in 

m 

in 

r- 

r- 

r- 

r- 

r- 

f" 

00 

OO 

00 

00 

00 

On 

On 

On 

On 

On 

00 

00 

OO 

00 

oq 

cq 

p 

p 

p 

p 

OH 

HH 

HH 

hh 

HH 

VO 

VO 

VO 

VO 

VO 

C4 

VO 

VO 

VO 

VO 

VO 

o 

o 

o 

o 

O 

hH 

O 

O 

O 

O 

O 

+ 

+ 

4" 

4- 

4" 

+ 

4- 

4- 

4" 

4- 

(2 

U 

UJ 

ta 

U 

PCI 

PCI 

PCI 

03 

03 

00 

OO 

OO 

00 

00 

CO 

CO 

CO 

CO 

CO 

OO 

00 

OO 

OO 

00 

in 

m 

in 

in 

m 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

cn 

ON 

m 

On 

CO 

ON 

CO 

ON 

CO 

On 

5 

5 

5 

5 

5 

— H 

HH 

O'! 

m 

Tf 

m 

— 

CM 

CO 

Tf 

m 

in 

m 

m 

in 

in 

O 

o 

o 

o 

o 

4* 

4- 

4- 

+ 

+ 

03 

03 

03 

03 

03 

00 

00 

00 

00 

00 

VO 

VO 

VO 

VO 

VO 

00 

00 

OO 

00 

00 

Tf 

Tf 

Tf 

Tf 

Tf 

♦ 

hh 

hH 

HH 

HH 

hh 

OO 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

m 

Tf 

o 

o 

o 

o 

o 

CM 

4- 

4* 

4* 

4- 

4- 

03 

03 

03 

03 

03 

CO 

CO 

CO 

CO 

CO 

O 

O 

O 

o 

O 

cd 

On 

ON 

ON 

On 

CfN 

m 

CO 

CO 

CO 

CO 

CO 

m 

\D 

p 

p 

p 

p 

CO 

in 

in 

in 

in 

iri 

in 

m 

m 

in 

in 

in 

o 

o 

o 

o 

o 

o 

4- 

4" 

4- 

4- 

4- 

4- 

03 

03 

03 

03 

03 

03 

OO 

00 

00 

OO 

OO 

OO 

CO 

CO 

CO 

CO 

CO 

OO 

On 

ON 

ON 

CTN 

On 

in 

CM 

CM 

CM 

CM 

CM 

CM 

hH 

hh 

hH 

hH 

hH 

CM 

VO 

vd 

vd 

vd 

vd 

H^ 

m 

m 

in 

m 

m 

in 

o 

o 

o 

o 

o 

o 

4" 

+ 

4- 

4- 

4* 

4- 

03 

03 

03 

03 

03 

03 

CO 

CO 

CO 

CO 

CO 

CO 

r- 

r- 

CM 

c^ 

c^ 

On 

On 

CTN 

VO 

VO 

VO 

VO 

p 

VO 

rq 

VO 

vd 

vd 

vd 

vd 

H^ 

m 

m 

m 

in 

m 

o 

in 

o 

o 

o 

o 

o 

o 

+ 

4* 

4- 

+ 

4* 

4- 

03 

03 

03 

03 

03 

03 

00 

00 

00 

OO 

OO 

OO 

8 

8 

8 

8 

8 

in 

VO 

hH 

hH 

hH 

hh 

On 

O 

hH 

hH 

hh 

hh 

hH 

CO 

CM 

C^ 

CM 

in 

m 

in 

in 

in 

in 

o 

o 

o 

o 

o 

o 

4* 

4" 

4- 

4- 

4* 

4" 

03 

03 

03 

03 

03 

00 

03 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

O 

5 

O 

VO 

O 

o 

o 

ON 

VO 

p 

p 

VO 

VO 

VO 

r- 

CM 

in 

in 

m 

m 

in 

CO 

m 

o 

o 

o 

o 

o 

o 

4- 

4- 

4* 

+ 

4" 

4- 

03 

03 

03 

03 

03 

03 

00 

OO 

OO 

OO 

OO 

OO 

r- 

r-- 

t-* 

r- 

CM 

o 

o 

o 

o 

o 

CO 

r- 

On 

On 

On 

ON 

2^ 

OO 

O 

p 

p 

p 

O 

.00 

00 

od 

od 

od 

cd 

in 

m 

in 

in 

m 

in 

o 

o 

o 

o 

o 

in 

o 

4- 

4- 

4- 

4* 

4* 

CO 

4* 

03 

03 

03 

03 

03 

03 

CO 

CO 

CO 

CO 

CO 

CO 

VO 

r- 

OO 

OO 

OO 

00 

OO 

r-* 

in 

p 

p 

p 

in 

Tf 

p 

od 

od 

od 

od 

CXD 

CO 

cd 

in 

m 

m 

m 

m 

m 

o 

o 

o 

o 

o 

o 

4" 

4- 

4- 

4- 

4- 

4* 

03 

03 

03 

03 

03 

03 

OO 

OO 

OO 

00 

OO 

CO 

OO 

Tf 

Tf 

Tf 

Tf 

CO 

On 

r- 

r- 

r- 

VO 

o 

o 

o 

o 

O 

hH 

On 

On 

CTN 

o^ 

ON 

Tf 

m 

m 

in 

in 

m 

CM 

in 

O 

o 

O 

o 

o 

CO 

o 

4* 

4- 

4- 

4- 

4- 

4- 

03 

03 

03 

03 

03 

03 

CO 

CO 

CO 

CO 

CO 

CO 

OO 

OO 

OO 

00 

00 

CO 

. 

00 

VO 

VO 

VO 

VO 

VO 

in 

m 

in 

m 

p 

p 

CO 

p 

On 

ON 

On 

On 

On 

Tf 

hh 

CM 

CO 

Tf 

m 

HH 

Appendix  A  SWIFTM86  Simulations 


A263 


Appendix  A  SWIFT/486  Simulations 


Appendix  A  SWIFT/486  Simulations 


A265 


VO  VO  VO  VO 

o  o  o  o 

W  W  W  W 
fn>r)COTt- 
^  O'  '  m 
^  oo  m 
^  oo 

rn  00  m  VO 
CO  Tt  cn  ^ 
VO  VO  VO  VO 

o  o  cp  o 
UJ  W  (JQ  t6 

VO  OO  VO  VO 
<N  Tt  VO 
CS  OO  <N  VO 
VO  OO  VO  «o 
rf  o  Tf 
lo  00  lo  cs 
O  »0  VO  VO 

9000 
1]  t6  uj  uj 

O  O  «0  VO 
D\  00  OV  VO 

'S-  o  Tt 
D  CN  O  O 
N  CS  CN  -H 
»  ^  00  Tf 

o  >o  10  VO 

9000 
ij  U  UJ  u 

O  <N  VO 
O  10  O  00 
D  Ov  O  »0 
»o  ^  ^ 
-j  vq  ^  VO 

^  VO 

n  VO  VO  VO 
D  O  O  O 

u  pb  uj  u 

o  Tt  cn  00 
^  0\  Ov  <N 
*0 

n  O  m  00 
n  o  m  VO 
-I  <s  VO 
VO  VO  VO 

p  o  o  o 
a  u  ub  u 

Ov  VO  CS 
}•  »-*  Tf 

r»  m  VO 
^  o 

f  ^  Tt  00 

-I  c4  ^  vd 

■)  VO  «0  VO 

p  o  o  o 
b  ub  ub  ub 

^  o\  00 
•)  <N  m  fo 
4  <N  <N  VO 

3  ^  VO  Tt 

4  0\  CS  »-« 

H*  ^  vd 
3  VO  VO  VO 

P  o  o  o 

^  ub  ub  ub 

D  00  00  Ov 

■>  VO  VO  ^ 

>  m  o  »o 

^  cn  to  VO 
*  rr  Tf  VO 
\  ^  Ov  Tt 

>  VO  VO  VO 

5  0  0  0 

I  pb  ub  pb 

I  rvj  ON 

>  00  CO  VO 
'  ^  0\  CO 

)  Ol  00  ov 
j  vq  »q  vq 
i  00  vd  ts 

I  VO  VO  VO 

>  cp  o  o 
I  pb  pb  pb 

'  Tf  TJ-  t3V 
I  Ov  VO 
tN  Tt  VO 

'  O  10 

VO 


«  CS  CO  Tf  VO 


00000 

u  pb  pb  pb  pb 

ts  ^  00  ^  cs 

Tf  VO  00  VO  rf 

VO  VO  CS  VO  VO 

^  0\  <N  Ov 

O  Tf  Tt  O 

CO  vd  vd  CO 

O  0\  Ov  O'  o 
^  o  o  o 

pb  ub  pb  u  pb 

Ov  ts  VO  CS  C7\ 
Qv  rf  VO  Tf  OV 
O  Ov  tS  On  O 
rfiOTf-voTt 

p  vq  cq  vq  o 

Ov  CS  ^  ON 
On  On  On  0\  0\ 

00000 

pb  pb  pb  pb  pb 

VO  00  Ov  00  VO 
O  00  rf  00  O 
VO  00  VO  00  VO 
CT\  On  O  ON  ON 
VO  00  VO 

vd  <s 
C^  00  00  00  ON 

00000 

pb  pb  pb  pb  pb 

Tf  CM  ^  CS  Tf 
On  VO  »0  »0  ON 
CO  VO  00  VO  CO 
<S  00  fS  ^ 
^  CO  00  CO  ^ 


)  00  00  00 
)  o  o  o 

]  pb  pb  pb 

>  ^  VO  00 

>  VO  ON 
)  O  VO  VO 
)  0\  00  VO 
;  p  Tt  00 

>  Tf  CO 

)  00  00 

>000 

I  pb  pb  pb 
■  On  00 
r-  r^i 
(  ^  CO  Q 
'  VO  O  O 
.  (N  NO 


J  tU  PJ  PQ 

3  »0  O 

F  00  ^ 

h  00  r- 

■>  VO  VO 
:>  On  O  O 

4  <N  ts  ^ 
'  r-  r-  r- 
p  o  o  o 

b  pb  pb  pb 

D  o  00  r- 

1  Tf  CO  00 

On  10  10 

5  <N  VO 

N  p  vq  CO 

F  vd  Tf  ts 

-  VO 

>000 

]  pb  pb  pb 

)  p  00  o 

^  O  CO  CO 

>  VO  00 

>  00  o  o 
)  cq  vq  ON 
V  On  rf 

>  VO  VO 

>000 

)  pb  pb  pb 

<  VO  rf 
'  On 
k  CO  CO 
’  O  Tt  ^ 

I  rq  00  ■'J 

I  r4  On 


pa  u  pj  pj  pj 
VO  00  CS  00  VO 
00  VO  VO  00 
ts  r4  <N  <N 
CO  VO  O  VO  CO 
Ttcsr^cSTj- 

K  ^  K 

VO  VO  VO  VO  VO 

pb  pb  pb  pb  pb 

VO  Tf  Tf  VO 
CO  On  fS  ON  CO 
CO  ON  ON  CO 
C7\  00  VO  00  On 

o  cq  ^  CN  p 

CO  vd  vd  CO 

pb  pb  pb  pb  pb 

VO  ^  NO 
00  VO  CO  VO  00 
CO  00  CM  00  to 
«0  VO  to  VO  «o 
cq  p  ^  cq 

— ^  CS  CN  CN  — ’ 
TfrfCOTfTf 

pb  pb  pb  pb  pb 

o  ^  o  o 

^TfONTj-rl- 

p  vq  rq  iq  ON 
rj-’  06  ^  00  rf 
to  CO  to  to  to 

pb  pb  pb  pb  pb 

VO  Tf  ON  rf  VO 

^  o 

^  VO  00  VO 
ON  O  CS  O  On 
p  CO  VO  to  00 
^  CO  to  ^ 
to  CN  CS  CN  CO 

pb  pb  pb  pb  pb 

O  CO  O  Tt 
>0  On  CO  ON  VO 
to  cs  to 

SCO  O  CO 

cq  r-  <s  <D 

r-:  ^  ^  ^ 

cs  CN  cs  cs  cs 

pb  pb  pb  pb  pb 

ON  CO  VO  CO  ON 
cs  ^ 

CO  ON  o  O'  CO 

vq  Tt  CN  Tt  VO 
CN  Tt  vd  Tf  CN 
cs  ^  —  cs 

pb  pb  pb  ub  pb 

00  CTn  CS  ON  00 
ON  Tt  -tt  Tf  C7\ 
CO  VO  VO  VO  CO 
Tf  00  ^ 

o<D  vq  10  00 

06  CN  0C5 


pa  ua  u  uj  ua 

^  <N  ^ 

O  —  CO  ^  O 

10  VO  VO 

Tt 

p  cq  Tt  cq  On 
cN  vd  vd  <s 
^  o  o  o  « 

pb  pb  pb  pb  pb 

VO  10  O  VO  VO 

'  CS  10  CN  ^ 

VO  Tl-  VO 

VO  ^  VO 


*—  CS  CO  ^  VO 


pa  u  pa  pa  pa 

VO  CO  to  VO 

CN  00  r-  00  CN 

CO  CS  to 

CN  r-  VO  tN 

Tf  cq  p  cq  Tf 

CN  CN  CS 
<N  ^  ^  CS 

CS  CN  CS  CS  <N 

pb  pb  pb  pb  pb 

VO  00  VO  ^ 
On  On  00  Ov  On 
ON  VO  CS  VO  ON 
r-  CO  00  CO 
P  ^  Tt  ^  O 

^  K 

CN  CS  CN  CS  <N 

pb  pb  pb  pb  pb 

to  CO  CN  CO  CO 
CN  00  O  00  CS 
VO  o  ^  «o 

00  On  00  ON  00 

vq  p  vq  p  VO 
CO  vd  vd  CO 
00000 

<N  CS  <N  CN  tS 

pb  pb  pb  pb  pb 

CO  CO  <N  CO  to 
CS  VO  O  VO  CS 
CN  VD  Tt  VO  C'i 
00  On  O  On  00 
p  cq  00  c» 

I-*  tN  CO  CS 


u  pa  pa  pa  tu 

-H  rt  VO  Tt  ^ 
Tt 

00  o  ts  o  00 
00  tN  tN  00 
vq  Tt  00  Tt  vq 

00  ^  cd 

On  Ov  On  Ov  On 

pb  pb  pb  pb  pb 

CO  00  00  CO 

CO  g  O  g  CO 

VO  CO  CO  CO  VO 
p  00  O  00  I-H 

Tt  vd  On  vd  Tt 
00  00  00  00  00 

pb  pb  pb  pb  pb 

CO  10  «o  10  CO 
Tt  00  O  00  Tt 
C'J  »0  C7N  VO  CN 
Tt  ^  VO  ^  Tt 
9  p  p  p  p 
^  cd  p  cd  •”< 
00  r-  r-  00 

pb  pb  pb  pb  pb 

VO  On  VO  On  10 
00  CO  00 

to  ^  Tt  ^  CO 

c?\  00  00  O' 

cq  Tt  p  Tt  00 
00  ^  ^  ^  00 
r-* 

pb  ub  pb  pb  pb 

00  ^  00 
VO  to  CO  VO 
^  VO  CJN  VO 

00  <s  00 
p  vq  p  VO  Ov 
cd  vd  oc)  vd  cd 
VO  VO  VO  VO  VO 

ub  pb  pb  pb  pa 

VO  VO  VO  VO  VO 
CO  00  Tt  00  CO 
00  VO  00  « 
CO  Tt  CO  Tt 
P  ON  ON  On 

ri  cd  CN  ^ 

—  tS  CO  Tt  VO 


A266 


Appendix  A  SWIFT/486  Simulations 


NORMAL  TERMINATION  (rnME=  95  :TIME  =  3.4214E+05) 
CPU  elapsed  time  =  21.26  seconds 


Appendix  A  SWIFT/486  Simulations 


A267 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  rtoortmo  burden  for  this  coiiection  of  information  is  estimated  to  average  1  hour  per  r^nse.  including  the  time  for  reviewing 
eiathwin^nd  maintaining  the  data  needed.  ar>d  completing  and  reviewing  the  collection  of  information.  Servd  comments  regarding  this  ^r^n  estimate  or  any  other  asp^  of  this 
collection  of  information  Jncluding  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services.  nr^nVn?  Jefferson 

OavisHighway  Suite  1204  Arlington.  VA  22202-4302.  and  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Project  (0704-0188).  Washington.  DC  20503. 


1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE 

April  1995  _ 


4.  TITLE  AND  SUBTITLE 

Evaluation  of  SWIFT/486  Model  with  Analytical  Solutions 


6.  AUTHORtS) 

Mansour  Zakikhani 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 


3.  REPORT  TYPE  AND  DATES  COVERED 

Final  report  _ 


5.  FUNDING  NUMBERS 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


U.S.  Army  Engineer  Waterways  Experiment  Station 
3909  Halls  Ferry  Road 
Vicksburg,  MS  39180-6199 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


Technical  Report 
EL-95- 17 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


U.S.  Army  Environmental  Center 
Aberdeen  Proving  Ground,  MD  21010-5401 


11.  SUPPLEMENTARY  NOTES 

Available  from  National  Technical  Information  Service,  5285  Port  Royal  Road,  Springfield,  VA  22161. 

12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 

12b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  is  unlimited . 

13.  ABSTRACT  (Maximum  200  words) 

This  report  describes  the  evaluation  of  SWIFT/486  by  comparing  computed  results  with  six  selected  analytical 
solutions  for  several  flow  and  solute  transport  scenarios  of  varying  boundary  conditions  and  solute  sources  in  porous 
media.  SWIFT/486  (Sandia  Waste-Isolation  Flow  and  Transport)  is  a  three-dimensional,  finite-difference  code 
which  can  be  used  to  simulate  steady-state  or  transient  flow  and  transport  of  chemicals  (including  brine  and 
radionuclide)  and  heat  in  porous  or  fractured  geologic  media.  The  geologic  media  may  be  homogeneous,  isotropic, 
heterogeneous,  and/or  anisotropic.  The  transport  processes  which  may  be  modeled  by  SWIFT/486  include 
advection,  dispersion,  sorption,  decay,  and  leaching.  Fluid  flow  of  variable  densities  and/or  viscosities  also  may 
be  modeled  by  SWIFT/486.  Either  a  radial  or  Cartesian  coordinate  system  can  be  used  for  domain  discretization. 
The  present  version  of  SWIFT/486  is  classified  as  a  single  phase  and  saturated  flow  model.  The  evaluation 
performed  here  complements  previous  SWIFT  evaluations  and  applications.  The  model  also  was  reviewed  for 
efficiency  of  coding,  convenience  of  input/output,  program  portability,  and  available  diagnostic  messages.  Note 
that  although  only  part  of  the  above  evaluation  steps  are  described  in  detail  in  this  report,  the  conclusions  for  all 
are  given.  Overall,  SWIFT/486  is  a  relatively  efficient  code,  requires  optimal  amount  of  computer  storage,  and 
has  sufficient  diagnostic  flags.  SWIFT/486  simulations  matched  closely  the  analytical  solutions  to  several  simplified 
problems. 


15.  NUMBER  OF  PAGES 

I  314 

Groundwater  Modeling  ■ 

Model  evaluation  SWIFT/486 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION  20.  LIMITATION  OF  ABSTRACT 
OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 


UNCLASSIFIED 


UNCLASSIFIED 


Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std.  239*18 
298-102 


